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Abstract. The stability of density-stratified viscous Taylor-Couette flows is considered using the Boussinesq 
approximation but without any use of the short-wave approximation. The flows which are unstable after the 
Rayleigh criterion (/t < f/ 2 , with fi — fi ou t/Oi n and f) — R ln /Rout) now develop overstable axisymmetric Taylor 
vortices. For the considered wide-gap container we find the nonaxisymmetric modes as the most unstable ones. 
The nonaxisymmetric modes are unstable also beyond the Rayleigh line. For such modes the instability condition 
seems simply to be (i < 1 as stressed by Yavneh, Mc Williams & Molemaker (2001). However, we never found 
unstable modes for too flat rotation laws fulfilling the condition fi > fj. The Reynolds numbers rapidly grow to 
very high values if this limit is approached (see Figs. and 0J. Also striking is that the marginal stability lines 
for the higher m do less and less enter the region beyond the Rayleigh line so that we might have to consider the 
stratorotational instability as a 'low-m instability'. 

The applicability of these results to the stability problem of accretion disks with their strong stratification and 
fast rotation is shortly discussed. 

Key words. Accretion, accretion disks, Hydrodynamics, Turbulence 
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1. Introduction 

The flow pattern between concentric rotating cylinders 
with a stable axial density stratification was firstly studied 
by Thorpe (1966) who concluded that stable stratification 
stabilizes the flow. The further experimental and theo- 
retical studies by Boubnov, Gledzer & Hopfinger (1995) 
and Caton, Janiaud & Hopfinger (2000) confirmed the 
stabilizing role of the density-stratification and showed 
that i) the critical Reynolds number depend on the buoy- 
ancy frequency (or Brunt- Vaisala, frequency) of the fluid 
and ii) the stratification reduces the vertical extension of 
the Taylor vortices. The computational results of Hua, Le 
Gentil & Orlandi (1997) have indeed reproduced the ex- 
periment results. 

The common feature of these studies is that the outer 
cylinder is at rest and flow is unstable after the Rayleigh 
condition for inviscid flow (which was extended to strati- 
fied fluids by Ooyama 1966), i.e. 



dR 



(R 4 n 2 ) < o 



(i) 



where Q is the angular velocity of the flow. 

Recently, Molemaker, McWilliams & Yavneh (2001) 
and Yavneh, McWilliams & Molemaker (2001) found 

AQ 2 

!m <0 < 2 > 



as the sufficient condition for (nonaxisymmetric) insta- 
bility. The condition (J2J) is identical with the condition 
for magnetorotational instability of Taylor-Couette flow 
(Velikhov 1959). These results have been derived by a lin- 
ear stability analysis for inviscid flow. The numerical re- 
sults of Yavneh, McWilliams & Molemaker (2001) demon- 
strate the existence of the hydrodynamic instability also 
for finite viscosity. 

An instability of density-stratified Taylor-Couette flow 
beyond the Rayleigh line ft = fj 2 for nonaxisymmetric dis- 
turbances have also been found experimentally (Withjack 
& Chen 1974). With their wide gap container (?) = 0.2) 
they found for fx > the stability curve crossing the clas- 
sical Rayleigh line. The observed instability was reported 
as nonaxisymmetric. The resulting experimental stability 
line, however, is very steep for positive fi (see their Fig. 8) 
and does never cross the line p, = fj. 

For real viscous flows there are very illustrative re- 
sults by Yavneh, McWilliams & Molemaker (2001). In 
the present paper a more comprehensive study of such 
flows is given. The governing equations and the restric- 
tions of the used Boussinesq approximation are discussed 
in Sect. 2 while the numerical results are presented in 
Sect. |21 Summary and final discussion are given in Sect. 01 

The existence of an instability in Taylor-Couette flows 
with a stable radial rotation law and with a stable z- 
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Fig. 1. The geometry of density-stratfied Taylor-Couette 
experiments. 



stratification of the density is a surprise in the light of the 
Solberg-H0iland criterion (see Riidiger, Arlt & Shalybkov 
2002). The necessary condition for stability reads 



dz 



dg z 
dR : 



(3) 



with g is the external force acceleration. Any conserva- 
tive force is a particular solution of J2J. Without external 
forces this relation is thus always fulfilled. If - as it is in ac- 
cretion disks - the gravity balances the pressure and the 
centrifugal force, then Eq. © is automatically fulfilled. 
Note that after the Poincare theorem for rotating media 
with potential force and Q = Q(R) both the density and 
the pressure can be written as functions of the general- 
ized potential so that © is always fulfilled. Generally, the 
magnetic field is not conservative and can never fulfill the 
condition This is the basic explanation for the exis- 
tence of the magnetorotational instability (MRI) driven 
by (weak) magnetic fields. 

Equation |J3J| has been derived by means of the short- 
wave approximation m < R/SR. For waves which are large 
in radial directions this condition might easily be fulfilled 
only for m = 0. It is, therefore, important to break the 
short-wave approximation in order to probe also the non- 
axisymmetric modes. For Kepler flows (with uniform grav- 
ity acceleration) this has recently been done by Dubrulle et 
al. (2004). With Boussinesq approximation and direct nu- 
merical simulation for viscous flow they find all stratified 
flows with negative &Q/&R unstable against nonaxisym- 
metric disturbances. According to their results a critical 
Froude number (as defined by Eq. I24fl exists below which 
the flow is stable. Whether the Boussinesq approximation 



can be used for too small Froude numbers seems still to 
be an open question. 

2. Equations and basic state 

In cylindric coordinates (i?, cf>, z) the equations of incom- 
pressible stratified fluid with uniform dynamic viscosity, 
p, are 
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+ (uv)iiji - 
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(5) 



(6) 



p is the density, P is the pressure, g is the gravity, v = p/p 
is the kinematic viscosity. 1 The equation which describes 
the evolution of the density fluctuation moving in the gen- 
eral density field is 



dp 
di 



{uV)p = 0. 



(7) 



We have to formulate the basic state with prescribed ve- 
locity profile u = (0, Rf2(R),0) and given density vertical 
stratification p — p(z). The system Q takes the form 



3. = ±.f*IL 

R pdR 1 
d 2 u rt> 1 du 4 



1 dP 

p dz 



= 0. 



dR 2 R dR R 2 
The last equation defines the angular velocity 
b 

R 2 ' 



O = a + — . 



(8) 



(9) 



where a and b are two constants related to the boundary 
values of the angular velocity, 0- ln , l? ou t, of the inner cylin- 
der with radius, i?i n , and the outer cylinder with radius, 
i?out- It follows 



a = 



p - r\- 



I 



b — n in R in 



rj 2 



1 



(10) 



the density diffusion term in the mass conservation equa- 
tion is neglected (see e.g. Caton, Janiaud & Hopfinger 2000) 
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with 
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(11) 



Differentiating the first equation of the system JHJ by z 
and the second equation by R, subtracting each other and 
using the supposed profiles of density and angular velocity 
one gets 



RQ 2 ^ 



0. 



(12) 



After this relation the density can depend only on the 
vertical coordinate z in the absence of rotation (Q = 0) 
and the angular velocity can only depend on radius in the 
absence of the vertical density stratification (dp/dz = 0). 
The supposed profiles of the angular velocity, ft = Q(R), 
and the density p = p(z) are, therefore, not self-consistent. 
Thus, we must admit more general profile for the density 
p = p(R, z) even though the initial stratification for the 
resting fluid is only vertical. In this case, the condition 
1121) takes the form 



2 dp ( ^dp 



0. 



(13) 



The fluid transforms under the centrifugal force from the 
pure vertical stratification at the initial state to mixed 
(vertical and radial) stratification under the influence of 
the rotation strongly complicating the problem. 

For real experiments the initial (without rotation) ver- 
tical stratification is small \d\ogp / dlogz <C 1 as is the ra- 
tio of centrifugal acceleration to the vertical gravitation 
acceleration 



R 2 n 



(14) 



so that after Ijl3|) the radial stratification is also small. Let 
us therefore consider the case of a weak stratification 



p +pi(R,z), pi<A), 



(15) 



where po is the uniform background density and l|13fl is 
fulfilled in zero-order. The perturbed state of the flow is 
described by 

ur, UA, + Rf2(R), u z , 

P {R) + Pi(R, z) + P, p + Pl (R,z) + p, (16) 

where |Pi/Po| <C 1 and ur, ua,, u z , P and p are the per- 
turbations. Linearizing the system (@J and selecting only 
the terms of the largest order we have 
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with uniform i/q = p/po- The first-order terms are left in 
the mass conservation equation due to vanishing of the 
zero-order terms. 

Due to (|14|) we can neglect dPo/dR in the first equa- 
tion and dpi/dR in the fourth arising from radial strat- 
ification (they will be \R 2 Q/g\ times smaller than terms 
arising from the vertical stratifications) and the system 
takes exactly the Boussinesq form 
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where N is the vertical buoyancy frequency with 



N z = - — 



9 dpi 
Po dz ' 



(18) 



(19) 



Suppose that the linear vertical density stratification 
dpi/dz = const and thus iV 2 is a constant, too. Then 
the coefficients of the system Ijl8(l only depend on the ra- 
dial coordinate and we can use a normal mode expansion 
of the solution F = F(i?)exp(i(m</> + kz + Lot)) where F 
represents any of the disturbed quantities. 

Let D — R out — Ri n be the gap between the cylinders. 
We use Rq = {RmD) 1 / 2 as the unit of length, the velocity 
Q m Ro as the unit of the perturbed velocity, Q m as the unit 
of w, N and Q. Using the same symbols for normalized 
quantities and redefining p as the dimensionless density 
pg/ poRo(2 2 n we finally find 
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dR 2 RdR V R 2 . 

— i Re(w + mfl)u z — i Re kP — Re p = 0, 
+ m(2)p- N 2 u z =0 (20) 

and 

flR ~R + * * ^ 
with the Reynolds number 

v 

The standard no-slip boundary conditions used at the 
inner and outer cylinder, i.e. 



Ur = US = u z = 0, 



(23) 



complete the classical eigenvalue problem. The same nu- 
merical method as in our previous papers about the 
Taylor-Couette problem (see e.g. Riidiger & Shalybkov 
2002) is used. Here we use a small negative imaginary part 
of u> to avoid problems with the corotation point to — mfl 
for m > 0. Thus, the calculated critical Reynolds numbers 
are not for the marginally stable state but for slightly un- 
stable state. To be sure that the calculated unstable state 
can be realized in experiments we checked the existence 
of the transition from stable to unstable state for several 
arbitrary points. 



imaginary part of oj is smaller than some fixed value de- 
pend on the vertical wave number. They have a minimum 
at a certain wave number for fixed other parameters. This 
minimum value is called the critical Reynolds number. 

In Fig. [5] we compare the calculated marginal stability 
line (i.e. ^s{uj) = 0) for axisymmctric disturbances with 
experimental values by Boubnov, Gledzer & Hopfinger 
(1995). The agreement is rather good except the small 
values of the Froude number 



Fr = 



f2[ n 
~N~ 



(24) 



This disagreement may indicate the violation of the 
Boussinesq approximation. For Kepler disks we find Fr ~ 
0.5 (Dubrulle et al. 2004). The unstratified fluids possess 
infinite Froude number. 



1000 



en 

m 800 



D 



if) 
Q 
_l 

O 

z 
>- 

UJ 

(Z 




600 
400 
200 b 



0.0 0.2 0.4 0.6 0.8 1.0 



50 c . . . i . . . i . . . i . . . i . . . ; 

O 2 4 6 8 1 O 

(N/O in ) 2 

Fig. 2. The marginal stability line for axisymmetric dis- 
turbances (m = 0) for f] — 0.78, fi — 0. The dots represent 
the experimental data of Boubnov, Gledzer & Hopfinger 
(1995). 

The code has been tested by computing for m — 1 the 
critical Reynolds number for the run 2 of Withjack & Chen 
(1974, their Tab. 1) with the experimental value 196.2 
(with our normalizations) and our computed result 200.6 
which we accepted to be in sufficiently good accordance. 

3. Results 

The imaginary parts of w, 3° (a;), decrease with increasing 
Reynolds number. The Reynolds numbers above which the 



Fig. 3. The marginal stability line for to = and critical 
Reynolds numbers for m > for fj = 0.78 and Fr = 0.5. 
The solid vertical line marks the p, = fj 2 limit, the dashed 
vertical line marks fl — f\ and the dotted vertical line 
marks fi = t) 1 ' 5 . 

The dependence of the critical Reynolds numbers on /t 
is given by Fig.|21for a narrow gap and in Fig.0]for a wide 
gap. The exact line of marginal stability is plotted only for 
to = 0. The axisymmctric disturbances are unstable only 
for p, < fj 2 in accordance to the Rayleigh condition (JTJ. 
For to > the slightly unstable lines with $s(uj) = — 10~ 3 
are given. 

The nonaxisymmetric disturbances are unstable also 
beyond the Rayleigh line (plotted as solid in the figures). 
The higher the to, however, the more the corresponding 
instability line approaches the Rayleigh line. Note, there- 
fore, that the 'stratorotational instability' (SRI, Dubrulle 
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Fig. 4. The same as Fig.|3but for a wide gap with fj = 0.3. 
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Fig. 5. The same as Fig. |3but for vertical wave number 
(left) and the pattern speed K(w)/to for to > (right). 



et al. 2004) only produces 1ow-to modes. This is an indi- 
cation that indeed within the short-wave approximation 
(high-m) it does not exist (see Riidiger, Arlt & Shalybkov 
2002). 

For nonstratified Taylor-Couette flows the nonaxisym- 
metric modes are only the most unstable disturbances 
for counter-rotating cylinders (e.g. Drazin & Reid 1981). 
The nonaxisymmetric instability of the stratified Taylor- 
Couette flow beyond the Rayleigh line (ft > fj 2 ) leads to 
the existence of some critical value, ft c , beyond which the 
nonaxisymmetric disturbances are the most unstable. Our 
results show that ft c ~ 0.27 and almost independent of the 
Froude number for fj = 0.78 (small gap) and /t c < for 
i) = 0.3 (wide gap). It is possible that for some fj the non- 



axisymmetric disturbances are the most unstable for all 
values of ft corresponding to unstable flows. 

There is another observation with the Figs. 01 and 
Approaching the line ft = fj the instability lines become 
more and more steep. We did not find any solution for ft > 
fj. If this is true one has to apply rather steep rotational 
profiles (not so steep as for nonstratified fluids but also not 
so weak as for the MRI) to find the modes of the linear 
SRI. This result seems to be of relevant for the discussion 
of the stability or instability of Kepler disks. 

For the narrow gap (fj — 0.78) the critical Reynolds 
numbers of the nonaxisymmetric modes only slightly de- 
pend on m. The same is true for the critical vertical wave 
number and the pattern speed !R(w)/to (Fig. - The ver- 
tical wave number only weakly depends on ft and the val- 
ues 5R(w)/to linearly run with ft. The situation is changed, 
however, for the wide gap (fj = 0.3). All parameters now 
strongly depend on both to and ft (Fig. ■ The trend for 
the vertical wave numbers is opposite for to = 3 to those 
for to = 1 and m = 2. 

The vertical wave numbers for both containers are of 
order 10 for to = 1. With our normalization the vertical 
extent of the Taylor vortices is given by 

" ^VW^F)- (25) 



R 



7T 

k 



With the mentioned value of k it is order of 0.1 for both 
the small-gap case and the wide-gap case. The cell be- 
comes thus rather flat. For nonstratified TC-flows one 
finds dz ~ i? ut — ^in while the cells under the influence 
of an axial magnetic field become more and more pro- 
late. The stratification generally reduces the height of the 
Taylor vortices. 
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Fig. 6. The same as Fig. [S]but for fj = 0.3. 

Unlike to the nonstratified Taylor-Couette flow the 
5ft(w) is not zero for stratified Taylor-Couette even for 
axisymmetric disturbances (to = 0). The onset of insta- 
bility is thus oscillatory ('overstability'). The question is 
whether a critical Froude number exists corresponding to 
the transition from stationary solutions to oscillating so- 
lutions? The answer is No. One cannot fulfill Eq. 120|) for 
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marginal stability (Q(u>) = 0) without a finite real part of 
u> for TV 2 ^ 0. The axially stratified Taylor-Couette flow 
bifurcates from the purely azimuthal flow through a direct 
Hopf bifurcation to a wavy regime. Depending on the value 
of p and fj this new regime can be either oscillating and 
axisymmetric or nonaxisymmetric and azimuthally drift- 
ing (see Figs. ISIEJ- For both our containers the pattern 
speeds are negative for positive m. The drift of the spirals 
is thus always in the direction of the cylinder rotation. 

Experiments have really demonstrated the oscillating 
onset of the axisymmetric instability (e.g. Caton, Janiaud 
& Hopfinger 2000). It would be interesting to design ex- 
periments with either rotating outer cylinder or wider gap 
to probe the bifurcation from the overstable oscillating 
axisymmetric flow pattern to the spiral nonaxisymmetric 
flow pattern. 

As an example for the container with the narrow gap in 
Fig. [7| the velocity eigenfunctions are presented for m = 1 
and for p = 0.7 exceeding the value of ff. The functions 
are smooth enough and do not suggest that the instability 
must be explained as a boundary effect. 
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Fig. 7. The velocity eigenfunctions for m = 1, fj = 0.78, 
p = 0.7, Fr = 0.5 with Reynolds number, vertical wave 
number and 3?(u;). The dotted lines are the real part and 
solid lines are the imaginary part. 
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Fig. 8. The Reynolds number, vertical wave number and 
5J(cj) in the transition from positive to negative for 
m=l,fj = 0.78, p = 0.7, Fr = 0.5. 

The Reynolds number, the vertical wave number and 
the real part of u> for the transition from positive to neg- 
ative imaginary part of u are given in Fig. GO i.e. for the 
transition from negative to positive growth rates. We find 
a continuous transition across the marginal stability line. 
It should thus be possible to realize the transition from 
stable to unstable flows in experiments. The vertical wave 
number and real part of uj are hardly influenced by the 
transition but, not surprisingly, the Reynolds numbers has 
a remarkably clear trend. 

4. Discussion 

It is shown that the Boussinesq approximation yields non- 
axisymmetric disturbances with low m of the stratified 
Taylor-Couette flow as unstable even beyond the Raylcigh 
line p, > rj 2 . Our results, however, also show that the criti- 
cal Reynolds numbers are extremely increasing approach- 
ing the line p, = r) so that as the condition for instability 
now the relation 



p < fj 



(26) 



seems to appear rather than p < 1 according to Yavneh, 
McWilliams & Molemaker (2001) and Molemaker, 
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McWilliams & Yavneh (2001). It is challenging to inter- 
prets the line p, = fj with the (galactic) rotation profile 
u$ =const in the same sense as to intcrprete the line 
fi = fj 2 with the rotation law for uniform specific angu- 
lar momentum R 2 Q = const. Below we shall consider the 
line p, = fj 15 as concerning the Kepler flow. 2 Dubrulle et 
al. (2004) for their model of rotating plane Couette flow 
have found unstable solutions also beyond the line jx = fj 
but also in these computations the rotation profile must 
be steeper than i? -2 / 3 (their Fig. 7). 

The SRI also leads to the situation that nonaxisym- 
metric disturbances can be the most unstable modes not 
only for counterrotating cylinders but also for corotat- 
ing cylinders. The characteristic values of fi where the 
nonaxisymmetric disturbances are the most unstable ones 
strongly depend on the gap width. For f) = 0.3 all positive 
jl are concerned (see Fig. 0J. It cannot be excluded that 
the nonaxisymmetric disturbances are the most unstable 
ones for all jl for f) smaller than some critical value. 

These results were obtained with the Boussinesq ap- 
proximation so that two restrictions remain. The vertical 
density stratification should be weak enough and the ro- 
tation should be so slow that the centrifugal acceleration 
can be neglected in total. If one of these conditions is vi- 
olated the Boussinesq approximation cannot be used and 
the situation becomes much more complicated. The dis- 
agreement between the calculated and the observed criti- 
cal Reynolds numbers for small Fr (see Fig.[5J may already 
indicate the violation of the Boussinesq approximation for 
strong stratifications. 

We have shown that for not too small negative d/2/di? 
Taylor-Couette flows with vertical density stratification 
become unstable against nonaxisymmetric disturbance 
with to = 1 even if they are stable without density strat- 
ification. Kepler flows seem to be concerned by this phe- 
nomenon. In the Figs. El E] the dotted lines represent the 
limit fi = fj 1 - 5 which might mimic the radial shear in 
Kepler disks. In both Figures we find a critical Reynolds 
number of only about 500 for the lowest (to = 1) mode. 
This is indeed a rather small number whose meaning, how- 
ever, should not be overestimated. Approaching the line 
fi = fj also these values become more and more large. 
Even more important is the finding that in magnetohydro- 
dynamic Taylor-Couette experiments the MRI only needs 
magnetic Reynolds numbers of O(10). This leeds to hy- 
drodynamic Reynolds numbers exceeding O(10 6 ) only for 
experiments with liquid metals in terrestrial laboratories. 
For hot plasma with magnetic Prandtl numbers of order 
10 (Noguchi & Pariev 2003) the necessary hydrodynamic 
Reynolds number is also only 0(1) or even smaller! 

The existence of the SRI might be important for as- 
trophysical applications. As suggested first by Richard & 
Zahn (1999) one should not forget (in particular for proto- 
planetary disks) to probe hydrodynamical instabilities as 
the source for the necessary turbulence in accretion disks. 

2 note, however, the general difference of the Kepler rotation 
law Q ex R~ 15 and the Taylor-Couette rotation law © 



In partial confirmation of results of Dubrulle et al. (2004) 
we have here demonstrated that for the rotation laws of 
Taylor-Couette flows under the presence of vertical den- 
sity gradients indeed linear hydrodynamic instabilities for 
low m exist. The properties of these modes are described 
above. Whether they are important for the accretion disk 
physics is still an open question. Note that the density 
stratification in accretion disks completely vanishes in the 
equatorial region and that the unstable modes discussed 
above would more or less lead to a nonaxisymmetric struc- 
ture of the disk rather than to turbulence. The next step 
in this direction must include the numerical simulation for 
well-designed but simplified global accretion disk models. 
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